Decoherence and the Nature of System-Environment Correlations 
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We investigate system-environment correlations based on the exact dynamics of a qubit and its 
environment in the framework of pure decoherence (phase damping). We focus on the relation of 
decoherence and the build-up of system-reservoir entanglement for an arbitrary (possibly mixed) 
initial qubit state. In the commonly employed regime where the qubit dynamics can be described 
by a Markov master equation of Lindblad type, we find that for almost all qubit initial states inside 
the Bloch sphere, decoherence is complete while the total state is still separable - no entanglement is 
involved. In general, both "separable" and "entangling" decoherence occurs, depending on temper- 
ature and initial qubit state. Moreover, we find situations where classical and quantum correlations 
periodically alternate as a function of time in the regime of low temperatures. 

PACS numbers: 03.65.Yz,42.50.Lc,03.65.Ta 



I. INTRODUCTION 

The increasing possibilities and sophisticated ways to 
control larger and larger quantum systems require a thor- 
ough understanding of the dynamics of open quantum 
systems [1, 2]. As prime example, the phenomenon of 
decoherence is an effect of the environmental impact on 
quantum behavior of the system, leading to a - some- 
times rapid - decay of its coherences in a certain ba- 
sis [3-7]. 

With respect to quantum information processing, de- 
coherence is the main obstacle that needs to be overcome. 
However, the influence of the environment on an open 
quantum system is not always of disturbing nature. For 
example, an environment can also mediate quantum cor- 
relations between two qubits [8, 9], which implies that 
in the process system and environment must get quan- 
tum mechanically correlated. It is an interesting question 
under which circumstances and how exactly those corre- 
lations develop. Many interesting results about correla- 
tions and entanglement in quantum many body systems 
have been obtained in the last few years [10]. The issue 
of system-environment correlations is also of importance 
for the quantum foundations of thermodynamics [11-13]. 

Although there are methods to detect the formation of 
system-environment correlations at the level of reduced 
dynamics [14], a detailed analysis of the role of quan- 
tum correlations in decoherence phenomena lies beyond 
such an approach. To gain deeper understanding one 
surely has to investigate the total (sometimes thermal) 
state of the composite quantum system, which in fact 
has been done for a couple of model systems [12, 15- 
19]. Some of the authors find decoherence accompanied 
by system-reservoir entanglement, some of them report 
of certain conditions under which decoherence appears 
without. Still, we lack a complete picture of entan- 
glement dynamics including finite temperature environ- 
ments and arbitrary system initial states. 



Here we present an approach to this problem, based 
on a partial P-representation of the total state, investi- 
gating a purely decohering qubit. It enables us to make 
detailed, state and temperature dependent statements on 
the evolution of the correlations between system and en- 
vironment, which we relate to the decoherence process. 



II. OPEN QUANTUM SYSTEM MODEL FOR 
DECOHERENCE AND REDUCED DYNAMICS 



Describing pure decoherence (phase damping), we will 
assume a standard model [20, 21] H to t = H SYS + Hj 
H env , by coupling a qubit H sys - 

Ox ® (Ea=i 9\A + 1 

harmonic oscillators H cnv = X)a=i ru^xa\ax 
Here Q denotes the transition frequency of the qubit and 
the coefficients g\ describe the coupling strength between 
the qubit and each environmental mode of frequency uj\. 
Assuming the environment to be initially in the thermal 
state ptherm = Gicp[- H cnv / k B T] /Tv[exp[- H cnv / k B T]] at 
temperature T, and the total initial state to be the prod- 
uct ptot(O) = Psys(O) ® Pthorm, the following exact time 
local master equation for the reduced density operator 
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T>(t) = cxp[— iQt — J ^(s)ds]. Here pij repre- 
sents the initial state of the qubit, with Bloch vector 
r = (x,y,z) = (cr). For the following, it is essential 
to consider arbitrary (mixed) initial qubit states. This 
choice reflects experimental limitations in the prepara- 
tion process, and also allows us to investigate decoherence 
and qubit-environment correlations in cases where the 
qubit is part of an entangled many qubit state [26, 27]. 
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While pure initial qubit states will always lead to system- 
environment entanglement, even a slight decrease of ini- 
tial purity may have a significant effect on the nature of 
the ensuing system-environment correlations. 

In (1) and (2) we have introduced the time dependent 
decoherence rate j(t) which by means of the spectral den- 
sity of the environment J(w) = ^2x=o \9x\ 2 ^( u — w a) can 
be written as 



lit) 



<hjJ{Gj)coih.[hLu/2k B T]cos[u)s]. (3) 



For a time independent j(t) = 7, qubit coherences decay 
on the decoherence time scale Tdcc = 7 _1 , sometimes 
referred to as "T2" . For the general case considered here, 
we define the decoherence time Td cc through the relation 
J^ dec j(s)ds = 1. The off-diagonal elements of the qubit 
have then decayed to a fraction e _1 of the initially present 
quantum coherences. 

Indeed, a constant decoherence rate j(t) = 7 is ob- 
tained for an Ohmic spectral density with a cut-off fre- 
quency uj c , in the limit of high temperatures ksT 3> fiuj c 
and large times uj c t ^> 1. Then eq. (1) is of Lindblad 
type, and the corresponding dynamics is Markovian. For 
our studies we choose a sharp cut-off, J(uj) = kujQ(oj — 
uj c ) [28], where k parametrises the coupling strength be- 
tween system and environment. Accordingly, the asymp- 
totic decoherence rate is given by — ATrnksT. 

Decoherence of qubits is sometimes modelled by ran- 
dom unitary dynamics [26, 29]. It is known that on 
the level of the reduced state, single-qubit decoherence 
can always be modelled that way, while for two qubits 
or more, there is genuine "quantum" decoherence - i.e. 
there is decoherence that - even on the reduced level - 
can only be modelled by coupling the system to a quan- 
tum environment [30-32]. Considering a single qubit 
and its given quantum environment, this raises the ques- 
tion to what extent the total state involves any quantum 
correlations. 

In fact, the most general class of total initial states 
which allows for completely positive reduced dynamics is 
classically correlated, with zero discord [33, 34]. We em- 
phasize that for our special choice of the spectral density, 
the map p ro d(i«) — > Pied{tf) induced by (1) is completely 
positive for all < U < tf. This is due to the positivity 
of 7(i) in our case, which also implies the " Markovian- 
ity" of the reduced dynamics for all times in the sense of 
both works [35, 36]. 

For all these recent investigations it is relevant to gain 
a deeper understanding of the build-up of quantum cor- 
relations between system and environment. In this paper 
we investigate the correlation dynamics of the total state 
for a mixed initial product state and identify entangled 
and classically correlated regimes. It turns out that in 
general the dynamics of system-environment correlations 
can be markedly rich. Let us provide the framework nec- 
essary to extract the desired information from the total 
state. 



III. THE TOTAL STATE 

We are interested in an exact expression for the to- 
tal state ptot(t), which we will represent in a coher- 
ent state basis. In a sense, our approach here is a 
mixed-state generalization of the non-Markovian quan- 
tum state diffusion approach to open quantum system 
dynamics [37, 38]. We will combine the environmental 
coherent state labels into the vector £ = (£1,^2, ■ ■ ■ ) of 
complex numbers and consistently make use of the nota- 
tion d 2 S,/n := d 2 £i/7r gP£ 2 /tt • • • (see also [38]). Further- 
more, introducing the average thermal occupation num- 
ber n\ of the A-th environmental mode, it is possible to 
expand the total state in terms of coherent states |£) 



Ptot(i) 



j^ e -^ /n p {t .^ n ^ ml (4) 



defining the matrix-valued partial P-function P(t) with 
values in the 2x2 dimensional state space of the 
qubit. Here we symbolically write cxp[— \t;\ 2 /n]/n := 

n A exp[-|6J7n A ]/n A . 

We find that expression (4) is a solution of the total 

von- Neumann equation with initial ptot(0) = /O sy s(0) <8> 

Ptherm if 

M+(t;£,r)Ax> B(t;£,OPoi 

where A ± = exp[-A(t) ± {(a(i)|f) + (f|a(t))}] and B = 
exp[-iSlt] exp[B(t) - {(b(t)\£) - (£[&(*))}]. Here we have 
introduced the complex time dependent vectors a(t) = 
(ai(i), a,2(t), • • • ) and b(t) with scalar product (a(i)|£) = 
J2\ a *\(t)£,\ an d vector components 

rt 
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Furthermore, we use the abbreviations 



A(t) =2Re[ ds [ dr V ^-\g x 
Jo Jo , n \ 



|2 g -iw A (t- 



B(t) = 2 Re / ds dr 
Jo Jo 



nx 



Initially, P = /5 sy s(0) = p rc d(0) and note that there are 
no approximations necessary to achieve the result (5) and 
thus via (4) to obtain the exact state of the composite 
system. 



IV. SYSTEM-ENVIRONMENT SEPARABILITY 
AND DECOHERENCE 

Having the total state at hand we now want to draw 
conclusions on its separability. For this purpose we argue 
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temperature [fico /ke] 



Figure 1. Separability of system and environment vs decoher- 
ence time - Lindblad regime: for temperatures ksT fko c 
the qubit may decohere without becoming entangled to its 
environment despite an initial purity of r = 0.98. 



as follows: 

As long as the partial P-function is positive-semi- 
dcfinitc the total state ptot(^) in representation (4) is 
trivially separable. Since one eigenvalue of P is always 
positive, a good indicator for this circumstance is the 
determinant dct[P] which in fact is independent of the 
coherent state labels £ and £*: 



det[P(t)] = e- 2A(t) [1 - z 2 ] 



,2B(t) 



(x 2 +y 2 ). (6) 



Recall that (x, y, z) are the Bloch coordinates of the ini- 
tial state of the qubit. This expression is positive as long 
as the following condition is fulfilled 



S(t) := A(t) + B{t) < i ln 



1-z 2 



x 2 + y 2 



(7) 



The quantity on the l.h.s. depends on the temperature 
of the surrounding heat bath via the thermal occupation 
numbers n\ = (ex.p[Huj / ksT] — and reads 



ft />S />OG 

S(T,t) = 4 ds dr duj 
Ja Ja Jo 



(8) 



x J(co) exp [hw/kT] cos[w(s — r)]. 



With our special choice of J(oj), S(T,t) can be written 
in terms of known special functions. For a given initial 
qubit state with (a) = (x,y,z), condition (7) defines an 
area in the (T, t)-plane, shown in fig. 1, within which the 
total state is necessarily separable. Note that for a pure 
initial state the r.h.s. of (7) is zero and the separable 
area vanishes. 

It turns out that as t — > oo, S(T,t) grows above all 
bounds for all temperatures and it becomes independent 
of temperature in the high temperature limit, S(T, t) — > 
S(t) w ln t + const. As a consequence, for any initial 




temperature [ftco c /kB] 



Figure 2. Fraction of qubit initial states that decohere while 
the total state is still separable: remarkably, in the regime of 
validity of the Lindblad equation (feT 3> fiu) c ) virtually all 
initial states decohere before any entanglement to the bath is 
built up. 



state and temperature, expression (4) of the total state at 
some point in time ceases to represent a proper mixture 
and thus, separability can no longer be concluded. Let 
us stress that positivity of P(t) is a sufficient but not a 
necessary condition for separability. 



Separability in the Lindblad Regime 

Condition (7) is visualized in fig. 1, where we relate our 
criterion for separability to the decoherence time Tdoc in- 
troduced earlier. We resort to the weak coupling regime 
(k = 10~ 3 ) and choose a qubit initial state in the equato- 
rial plane (z — 0) of the Bloch sphere with purity r — 0.98 
(recall that r denotes the length of the Bloch vector, i.e., 
r = |r|). 

While at low temperatures Td ec lies outside of our sep- 
arability area, it is inside this area for temperatures of 
about ksT > 2 fko c . In the high temperature regime, 
when the dynamics is governed by a Lindblad equation, 
decoherence of the qubit is thus complete, while the to- 
tal state is still separable. This is especially striking be- 
cause the initial state is so close to being the pure state 
(|0> + |l»A/2. 

The question arises as to how many initial states, de- 
pending on temperature, meet the very same fate. In 
fig. 2 we display the relative volume of the set of states 
that are still separable at Tdcc, i.e. after decoherence. Re- 
markably, in the regime where the master equation (1) 
is of Lindblad form (7 = const) for practically all states 
decoherence has long happened before any entanglement 
between system an environment can build up. 

In fact, this is not too surprising: in the weak coupling 
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limit ( "Born approximation" ) a master equation of Lind- 
blad type is usually derived assuming that the joint state 
of system plus environment to the lowest order in the in- 
teraction remains a product state [2, 39]: one substitutes 
Ptot(i) ~ Psys(t) ® Pthcrm in the perturbative evolution 
equation. The severity of this approximation is appar- 
ently manifested in the fact that entanglement builds up 
very slowly compared to the time scale given by decohcr- 
ence. 



V. SYSTEM-ENVIRONMENT 
ENTANGLEMENT AND DECOHERENCE 

All pure qubit initial states different from |0) and |1) 
will lead to system-environment entanglement immedi- 
ately. However, as we have shown in the previous section 
even a slight reduction of purity (e.g. r — 0.98) may 
lead to a prolonged period of separability, long enough 
for the qubit to loose its coherence. Being separable at 
least up to the time when det[P] becomes negative, p t ot 
turns into an entangled state at the latest when its par- 
tial transpose p\^ t yields a negative expectation value 
£ PT = (^Ip^J I*) in some state l^) of the composite 
system [40]. Hence, to detect entanglement we aim to 
find a state \ip(t;^*)) ~ (£ \^(t)} in the Hilbert space of 
the qubit such that 



7r n 

Denoting the pure state on the Bloch sphere 
which has the smallest overlap with the trans- 
pose of the initial state of the qubit as (u,v) := 
(--v/(l - z/r)/2,e^V(l + z/r)/2), it turns out that 



KHt;0> = 



e -K|o+6)/2+»nt/2 u 
_ e (£|a+6)/2-»fit/2 v 



(10) 



fulfills this task. This choice of state is further elaborated 
upon in section VI. By solving the integral (9) one finds 
the following expression for the expectation value: 



£^{t) = 



2(n+l) 



(,,-«>+¥ 


V--\ 




r 



?(t) {x 2 + y 2 ) 



(11) 



where S arises from S = A + B by inverting the coeffi- 
cients (n\ + l)/n\ — > n\/ (n\ + l). This quantity becomes 
negative, i.e. the total state is entangled, if the following 
condition is fulfilled: 



E(T,t) := S(T,t) - S{T,t) > In 



x 2 + y 2 



(12) 



Note that (r — z 2 ) / (x 2 +y 2 ) > 1 and hence the quantity 
on the r.h.s. (12) is always non-negative. 



As before with inequality (7) for separability, for any 
initial state, this new inequality (12) defines an area in 
the temperature-time plane in which entanglement of the 
total state is assured. As before, this is a sufficient but 
not a necessary condition for entanglement. 

In terms of the spectral density the relevant function 
reads 



/>t />S POO 

E(T, t) =8 ds dr div 
Jo Jo Jo 

x J((jj) sinh [fruj/kT] cos[w(s — r)]. 



(13) 



In the limit of high temperatures and large times, 
E(T, t) becomes independent of time, E(T,t) ->• E(T) w 
^Khwc/ksT. As a consequence, for every initial state 
there is a temperature above which our criterion is no 
longer able to detect entanglement. As we will show in 
figures 3 (a) and (b) for low temperatures and relatively 
short times the criteria for separability and entanglement 
cover almost the entire (T, t)-plane and only little "un- 
known territory" (the white areas) remains. 



Separability and Entanglement in the Strong 
Coupling Regime 

In fig. 3 we relate the decoherence time Td cc (T) to our 
criteria for separability and entanglement. To ensure that 
entanglement builds up fast, we here resort to the strong 
coupling regime (k — 1). It becomes apparent that for an 
initial qubit state that is pretty pure (r = 0.9), in the con- 
sidered low temperature range, decoherence is complete 
only after entanglement has built up (a). By contrast, 
for an initial qubit mixed state with purity r = 0.1 (de- 
pending on temperature) decoherence can occur before 
or after the two subsystems were able to entangle (b). 

These results show that in general there is no direct 
connection between decoherence and the formation of 
system-reservoir entanglement: surely, if the qubit is pre- 
pared in a pure state, decoherence is always accompanied 
by entanglement of the total state, independently of tem- 
perature. If the initial state is already a completely in- 
coherent mixture, the two subsystems will never get en- 
tangled. In all other cases there are two disjoint domains 
of temperature implying decoherence before and after the 
qubit was able to entangle with its environment, respec- 
tively. 

We now want to further elucidate the nature of deco- 
herence depending on the initial state of the qubit. In the 
cut through the Bloch sphere of fig. 4 all states inside the 
area enclosed by a curve of given temperature decohere 
while the total state is still separable. In similar vein, in 
fig. 5 all states outside the area determined by a curve of 
given temperature decohere after entanglement has built 
up. 

As the coupling to the environment distinguishes the 
south and north pole of the Bloch sphere as robust states, 
for a given purity r, there is no global value of temper- 
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(a) 




0.5 1 1.5 2 

temperature [fi© /k B ] 

(b) 




0.5 1 1.5 2 

temperature [fico c /kB] 

Figure 3. System-environment correlations vs decoherence 
time - strong interaction regime: (a) For initial states of high 
purity (here r — 0.9) decoherence mainly takes place after the 
total state has become entangled, (b) for mixed initial condi- 
tions decoherence can occur in both domains, here shown for 
an initial state with purity r = 0.1. 

aturc T defining the cut between "separable" and "en- 
tangling" decoherence. Rather, states close to the poles 
tend to decohere while the total state is still separable. 
States near the equator, by contrast, tend to decohere 
involving entanglement. 

Oscillations between Separability and Entanglement 

As already mentioned in the beginning, the dynamics 
of system-environment correlations can be quite rich. In 
the low temperature regime, in particular, we are able 
to detect conditions where the total state alternates be- 
tween being entangled and separable. In fig. 6 we show 
the correlation dynamics for an initial qubit state with 
purity r = 0.95. It becomes apparent that for tempera- 
tures around T = 0.3 hw c /kB the nature of the correla- 
tions changes oscillatorily between classical and quantum 




1 0.5 0.5 1 



Figure 4. Separable decoherence: cut through the Bloch 
sphere at (j> = 0. All states inside the area enclosed by a 
curve of given temperature decohere while the total state is 
still separable. Temperatures in units of &j c /fcs. 




Figure 5. Entangling decoherence: cut through the Bloch 
sphere at <p — 0. All states outside the area determined by a 
curve of given temperature decohere after entanglement has 
built up. Temperatures in units of hwc/ks- 

as a function of time. Although the special form of these 
oscillations depends on the form of the spectral density 
J(u>), we found numerically that similar behavior exists 
for different choices of J. 



VI. QUALITY OF ENTANGLEMENT 
DETECTION FOR A SINGLE MODE 
"ENVIRONMENT" 



Fig. 6 demonstrates that - at low temperatures at least 
- our separability (7) and entanglement (12) criteria are 
quite good: the "unknown territory" of correlations is 
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0.1 0.2 0.3 0.4 0.5 
temperature [fico c /kB] 

Figure 6. Oscillatory behavior in correlation dynamics: at 
low temperatures one can find conditions where separable and 
entangled domains alternate, here shown in the weak interac- 
tion regime (k = lCT 3 ) for a qubit initial state with purity 
r = 0.95 (z = cf> = 0). 




temperature [fico/k B ] 

Figure 7. Single mode environment - comparison to numerics: 
temperature dependence of the exact instant r cr i t (T) at which 
the total state changes from PPT to NPT, and vice versa 
(points). The result obtained by our criterion (12) is displayed 
by the shaded areas. Obviously, there is perfect agreement. 
Qubit initial state with purity r = 0.75, z-coordinate z — 0.2 
and azimuth angle <j) — 0. 

very small. To also estimate its accuracy for higher tem- 
peratures we study entanglement dynamics of the total 
state in the case of a single mode "environment", where 
the full dynamics can be accessed numerically. We are 
interested in the exact instant r cr it(T) at which the defi- 
niteness of the partial transpose pfj^ ceases to be positive 
(PPT) and negative eigenvalues emerge (NPT). On the 
other hand eq. (12) allows one to analytically derive the 



instant r ont (T) from which on we are able to detect en- 
tanglement with our test state determined from eq. (10). 

In fig. 7 we compare the temperature dependence of 
these characteristic times for an initial state with purity 
r = 0.75. It turns out that there is an excellent agree- 
ment between the exact numerical and the analytical re- 
sult, even at high temperatures. We are confident that 
also in the case of larger environments, our entanglement 
criterion (12) is able to very efficiently detect a negative 
partial transpose. Clearly, it is an interesting and open 
question, whether the white areas in our (T, i)-diagrams 
correspond to entangled, bound entangled or separable 
total states. 



VII. CONCLUSIONS 

We have investigated the dynamics and nature of 
system-environment correlations of a decohering qubit 
with arbitrary (mixed) initial state coupled to a bath 
of quantum harmonic oscillators. When the reduced dy- 
namics can be described by a Markov master equation of 
Lindblad type, we found that for almost all qubit initial 
states inside the Bloch sphere decoherence is complete 
long before the build-up of any system-reservoir entan- 
glement: the total state is still separable. 

By contrast, in the strong coupling regime the correla- 
tion dynamics was found to be more complex: dependent 
on initial state and temperature of the surrounding heat 
bath, decoherence can occur before or after the qubit is 
able to entangle with its environment. At very low tem- 
peratures we were even able to detect conditions where 
the total state alternates between separable and entan- 
gled domains as a function of time. 

Finally, we have compared our criterion for system- 
reservoir entanglement with exact numerical results for 
the occurrence of a negative partial transpose in the 
case of a single mode environment and found excellent 
agreement. We therefore believe to have found a reli- 
able criterion for system-bath entanglement. Clearly, it 
is desirable to further investigate the nature of system- 
environment correlations in the white regions in order to 
cover the whole range of temperatures and times. More- 
over, we are confident that our approach may be applied 
to more general open quantum system dynamics. 
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